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ABSTRACT 

Context. Extrapolations of solar photospheric vector magnetograms into three-dimensional magnetic fields in the chromosphere and 
corona are usually done under the assumption that the fields are force-free. This condition is violated in the photosphere itself and a 
thin layer in the lower atmosphere above. The field calculations can be improved by preprocessing the photospheric magnetograms. 
The intention here is to remove a non-force-free component from the data. 

Aims. We compare two preprocessing methods presently in use, namely the methods of Wiegelmann et al. (2006) and Fuhrmann et 
al. (2007). 

Methods. The two preprocessing methods were applied to a vector magnetogram of the recently observed active region NOAA 
AR 10953. We examine the changes in the magnetogram effected by the two preprocessing algorithms. Furthermore, the original 
magnetogram and the two preprocessed magnetograms were each used as input data for nonlinear force-free field extrapolations by 
means of two different methods, and we analyze the resulting fields. 

Results. Both preprocessing methods managed to significantly decrease the magnetic forces and magnetic torques that act through the 
magnetogram area and that can cause incompatibilities with the assumption of force-freeness in the solution domain. The force and 
torque decrease is stronger for the Fuhrmann et al. method. Both methods also reduced the amount of small-scale irregularities in the 
observed photospheric field, which can sharply worsen the quality of the solutions. For the chosen parameter set, the Wiegelmann et al. 
method led to greater changes in strong-field areas, leaving weak-field areas mostly unchanged, and thus providing an approximation 
of the magnetic field vector in the chromosphere, while the Fuhrmann et al. method weakly changed the whole magnetogram, thereby 
better preserving patterns present in the original magnetogram. Both preprocessing methods raised the magnetic energy content of 
the extrapolated fields to values above the minimum energy, corresponding to the potential field. Also, the fields calculated from the 
preprocessed magnetograms fulfill the solenoidal condition better than those calculated without preprocessing. 

Key words. Sun: magnetic fields - Sun: atmosphere - magnetohydrodynamics (MHD) 



' 1. Introduction ized by the equations 

t-h ; 

[ In order to model and understand the physical mechanisms un- V x B — a(r) B , (1) 

t-H derlying the various activity phenomena that can be observed V ■ B = (2) 

in the solar atmosphere, like, for instance, spots, faculae, flares, 

. i-H \ and coronal mass ejections, as well as their mutual connections where a(r) denotes a scalar function of position r which, because 

^ ■ an interactions, the magnetic field vector throughout the atmo- of Eq. ©, is constant along the magnetic field lines. As an im- 

5^ \ sphere must be known. Unfortunately, reliable magnetic field plication of the dominance of the magne tic energy density over 

measurements are still restricted to the level of the photosphere, all other energy densities (cf. lGarvll2001l) . this approximation is, 

where the inverse Zeeman effect in Fraunhofer lines is observ- except for the times of explosive events, presumably valid from 

able. Progress is very slow here, due to fundamental difficulties the upper chromosphere up to coronal heights of ~ \Rq above 

in unambiguously deriving the magnetic field from polarimetric the photosphere. 

measurements in chromospheric or coronal spectral lines. Other The development of extrapolation methods began with cal- 

information sources, like non-burst radio emission or contrast culations of potential fields, which correspond to the case of 

pictures in selected spectral lines or continuum parts of the solar a = 0, either abo ve limited photospheric areas , in pa rtic- 

spectrum, provide only order-of-magnitude estimates or limited ular active regions dSchmidtl Il964t iTeuber et all 1 19771). or 



qualitative information. in spherical shells above the full p hotosphere dSchatten et all 



As an alternative to measurements in the chromosphere and 119691: lAltschuler & Newkirki 1 19691) . Extrapolation methods us- 

in the corona, the measured photospheric field may be extrap- ing linear force-free fields, which correspond to the case of 

olated theoretically into these higher layers of the atmosphere, a spatially constant non-vanishing a, foll owed as a next step. 

Such extrapolations are usually done under the assumption that Corres p onding solutio n s wer e given b y |N akagaw a & Raadu 

the magnetic field B is approximately force-free, i.e. character- (Il972l) . IChiu & Hilton! dl977l) . ISeehafetl d 19781) . lAliss andrakis 
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(fl98l- 



■Semel (1 19881): see also rev i ews a n d compariso ns by 
Seehaferl (1 1 9821) . LSeehafer & Staudd (Il983h . iGarvl (Il989h . and 
Sakuraildl989l 

Over the last three decades, extrapolation methods for non- 
linear force-free fields were also developed. Here a is al- 
lowed to vary spatially. Unlike the extrapolation methods for 
current-free and constant-ff force-free fields, which require 
only line-of-sight magnetograms, the non-constant-or force- 
free fields are calculated from photospheric vector magne- 
tograms. Presently great efforts are made to improve the 



methods used here. The y include (i) direct vertical i ntegra- 
tion schemes (|Wu et all 1 19851 Il990t ICuperman etall 1 1990b 



iDemoulin et al 



19921: Song et al l 12006). (ii) G rad-Rubin meth- 



ods dGrad & Rubinb 11958b ISakuraib Il98ll lAly & Seehaferl 



fl99l lAmari et al. 1 119971119991120061 Wheatlan 
bound a ry-integral methods (lYan & Sakura 



dl 20041 
I TT9971 



20061 



20001 



„ 12003b lYan & LiL I2006L iHe & Wangb 120061). Civl relax- 
ation methods dYang et all 119861: iMikic & McClvmontl 11994 



Valori et al 




iRoumeliotisl 119961: 
tion methods (IWheatland et al 
120071) . Discussions and quantitative 
of the nonlinear extrap o lation 



and (v) op t imiza- 
IWiegeimannl 12004 



and 



m 



comparisons 
found, e.g 



tests 
methods 
Alvl rfl989h ISakurail (fl989l . 
ISchriiver et ail d2006l). iMetcalf et al.1 (120081) ISchriiver et al.l 



McClvm ont et al.l d 199 



(l2008h . lWiegelmannl(l2008h . and lDeRosa et al.l(l2009l). 

Ap plied to semi-analyt ic al forc e-free test fields (lLow & Loul 

[l990l iTitov & DemouM [l999h . the nonli near extrapola- 
tion methods show encouragi n g results dS c hriiyer et all 
12006b IWiegelmann etafl l2006ab IValori et all I2007L l2010h . 
Nevertheless, most of these methods are still lacking a sound 
mathematical basis. The existence of solutions seems to have 
been proven merely for Grad-Rubin methods, and here only for 
suffic i ently small \ a\, i.e. close to the potential field (Bineaul 
1 1972b iKaiser et all l2000b iBoulmezaoud & Amarib l2000h . The 
Grad-Rubin methods do not exploit the full information content 
of photospheric vector magnetograms (the normal field compo- 
nent B„ and a are used as boundary data, but a is prescribed 
only on a part of the boundary where B„ is either positive or 
negative). 

By contrast, the fields calculated by direct vertical integra- 
tion always satisfy the conditions given by the photospheric vec- 
tor magnetograms, and are valid and uniquely determined so- 
lutions to the problem, whenever such solutions exist and are 
sufficiently smooth. The solutions are however unstable, or the 
problem is ill-posed, in the sense that small changes of the pho- 
tospheric boundary values can lead to large changes of the field 
above the photosphere. 

The boundary-integral, relaxation and optimization methods, 
on the other hand, impose conditions on all three components of 
the magnetic vector on the complete boundary of the solution 
domain, mostly a rectangular cuboid or a spherical shell above 
the photosphere, with some assumtion for the non-photospheric 
parts of the boundary. This way the instability of the solution to 
boundary-value perturbations is avoided, but inconsistencies in 
the boundary values lead to only partially force- and divergence- 
free reconstructions. Therefore, the prescribed boundary values 
cannot be arbitrary but have to satisfy consistency conditions. 
This consistency of the data with the condition of force-freeness 
cannot be expected in practice, however, because (i) the field is 
not force-free in the photosphere and in the lower chromosphere 
(possibly with the exception of regions with particularly strong 
fields, like sunspots) and (ii) compatibility with the condition 
of force-freeness, if it were given for a real photospheric field, 
would be destroyed by measurement errors, which are always 



present, in particular in the transverse field components (perpen- 
dicular to the line of sight of the observer), e.g. caused by sys- 
tematic uncertainties (typically much larger for the transverse 
than for the longitudinal field components), or by an inappropri- 
ate resolution of the 180° ambiguity. 

The compatibility of the boundary data with the condition 
of force-freeness can be improved by an appropriate preprocess- 
ing of the dat a. Algorithms for such a pre proc essing have been 
sugges ted by IWiegelmann et al.l d2006bl) and iFuhrmann et al.l 
(120071) . Their strategy is to modify the magnetographic data 
within certain margins such as to minimize the total magnetic 
force and the total magnetic torque on the volume considered. 
These integral quantities vanish if the magnetic field is force- 
free, and ca n be expressed as a function of the bo undary values 
of B alone dMolodenskiib Il969t IaR fl984l [l989t ILowL [1981 . 
The minimization process removes a non-force-free component 
from the data, so that the resulting boundary values can be more 
consistently extrapolated into a force-free field. 

Another, related objective of the preprocessing is to remove 
small-scale noise from the data. It has been observed that the 
convergence of the evolutionary or iterative sequences calculated 
using the relaxation or optimization methods can worsen sharply 
due such a noise, and the quality of the obtained solutions, that 
is, the degree to which they are force-free, is also worsened. This 
phenomenon is not well understood presently. 

One aim of smoothing the photospheric magnetogram is also 
to mimic the expansion of the solar magnetic field between pho- 
tosphere and chromosphere, where the field is assumed to be- 
come force-free. The idea is that force-free models are only valid 
above the chromosphere, but not suitable for modelling the high- 
beta plasma between photospher e and chromosphere; see also 
discussion in lMetcalf et al.l(l2008l) . In this sense the preprocessed 
field may be considered as an approximation to the chromo- 
spheric field, which is much smoot her and shows les s fine struc- 
tures than the photospheric field. iJing et all J2010h compared 
unprocessed and preprocessed photospheric field measurements 
from Hinode with chromospheric field measurements from Solis 
and found that the preprocessed field is indeed a reasonable ap- 
proximation of the chromospheric field. 

In this paper, we com pare the preprocess ing methods of 
IWiegelmann et al.l (l2006bl) . hereafter ppTW, and lFuhrmann et all 
2007ir hereafter ppMF. In particular, the two methods are ap- 



plied to a vector magnetogram of a recently observed active re- 
gion, NOAA Active Region (AR) 10953, which has been the 
target in a comparitative study of different extrapol ation methods 
for non linear force-free magnetic fields reported in lDeRosa et al.l 
(120091) . 

The remainder of the paper is organized as follows: In Sect. [2] 
we describe the two preprocessing methods. Then, in Sect. [3] we 
apply them to AR 10953: in Sect. l3.1l some observational details 
of AR 10953 and of the data used are given, and in Sects. 13.21 
and !3.3l the effects of the preprocessing on the magnetogram and 
on the extrapolated field, respectively, are analyzed. In Sect. |4] 
we draw conclusions and discuss our results. 



2. Preprocessing algorithms 

2.1. General strategy 

The general strategy of ppTW and ppMF is to modify the bound- 
ary data given by the photospheric vector magnetogram in order 
to decrease the total force and the total torque on the volume 
considered and to reduce the amount of small-scale noise in the 
data. This is done by minimizing a functional L of the boundary 
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values of B. L is the sum of a number of subfunctionals each 
of which measures one of the quantities to be made small, i.e., 
the total force, the total torque, and the amount of noise. The 
degree of deviation of the modified photospheric fields from the 
original one during the minimization is controlled either via an 
additional subfunctional of L that measures the deviation, or by 
allowing the modification of the data only within prescribed do- 
main borders. The different subfunctionals of L are weighted in 
order to control their relative importance for the minimization. 



2.1 .1 . Total magnetic force and total magnetic torque 
Let 

T=\ jxBAV 

Jv 

denote the total magnetic force and 



JV 



\ rx(jxB)AV 
Jv 



(3) 



(4) 



the total magnetic torque on the volume V, the solution domain. 
In Cartesian coordinates x, y, z, with the z axis pointing upward 
in the atmosphere, and with S denoting the magnetogram area 
in the plane z = 0, the condition T = is approximated by 



m f x = - J B X B Z Ax Ay = , 
fiofy = - J B y B z dxdy = 0, 
2/4, T z = J (B 2 X + B] - B 2 ) Ax Ay = , 
while the condition N — is approximately expressed by 
2fi N x = J y (b 2 + B 2 x - Bfj Ax Ay = , 
2fi N y = J x (-B 2 . - B 2 + B 2 ) Ax Ay = , 
fxo N z = j (yB x B z - xB y B z ) Ax Ay = 



(5) 
(6) 
(7) 

(8) 
(9) 
(10) 



( iMolodenskiil 11 9691: IXR [19891: iFuhrmann et al.1 12001 . For ex- 
pressing the conditions of vanishing force and torque on V ex- 
actly, the integrals on the right-hand sides of Eqs. <f5T> — (fTOb would 
have to be extended to integrals over the complete surface dV 
of V. The restriction to integrations over the photospheric mag- 
netogram area S is done under the assumption that all relevant 
magnetic flux is closed on the photosphere, and the field on the 
rest of the boundary is so weak that its contribution to the surface 
integrals for T and N is negligible. 

2.1.2. Smoothing 

There are several methods for removing small-scale irregulari- 
ties from a noisy vecto r magnetogram, hereafter called smooth- 
ing dPress et al.L Il989h . We restrict ourselves to the discussion 
of the two smoothing methods employed by the codes discussed 
below. 

The first method consists in applying the two-dimensional 
Laplacian, A„. = d 2 /dx 2 + d 2 /dy 2 , to the photospheric field, 
B{x, y,z — 0), and minimizing the integral of the square of the 
obtained function, [A xy B(x,y, z = 0)] 2 , over the magnetogram 



area. [A„B(x, y, z — 0)] measures the roughness of the pho- 
tospheric boundary data, and minimizing it corresponds to rec- 
ducing the deviations of the field values at given points from 
the mean values of the field in points around them, in accor- 
dance with the mean-value property of the harmonic functions. 
Since the second-order derivatives measure the curvature of the 
surfaces B x (x, y,z = 0), B y (x, y,z — 0), and B z (x, y,z — 0), de- 
creasing the modulus of A xy B(x,y,z - 0) to zero on S would 
remove all global or local maxima or minima of the three field 
components (except for those on the boundary line of the mag- 
netogram), including the physically significant ones, and the to- 
tal energy content of the field might be strongly reduced. This 
kind of smoothing may however well mimic the expansion of 
the magnetic field between photosphere and chromosphere. 

The second smoothin g method consid ered here is the 
windowed-median method ([Press et allll989h . Compared to the 
first method, essentially the normal, arithmetic mean is replaced 
by the median. That is, one calculates the median of the field val- 
ues in a small window around a given grid point and then min- 
imizes the difference between the field value at the grid point 
and the median of the window. This way one tries to avoid a too 
strong influence of large disturbances that occur only in a small 
number of points. Also, the structures in the measured fields may 
be better conserved even if the quantity minimized here is de- 
creased to zero. In general, however, decreasing the total force 
and the total torque on the one hand and smoothing on the other 
hand are competing objectives, so that a too strong smoothing 
would prevent a sufficient reduction of the total force and the 
total torque. 

2.2. Preprocessing to mimic the chromospheric field fppTW) 

IWiegelmann et al" | d2006bj) developed a preprocessing routine in 
which the ideas described above are implemented. In discretized 
form, the square of the total magnetic force, as approximated by 
Eqs. ©-(EK is given by 



U = 





2 




2 / 


Yj B * B z 


+ 




+ 


V p 




\ p > 


V 



^{b x + b 2 -b 2 ) 



\ p 



(11) 



and the square of the total magnetic torque, as approximated by 
Eqs. ©-Coll, by 



y^yB x B z -xB y B z 



(12) 



where the summation is over the points P of a photospheric 
grid (in the numbering of the different subfunctionals we follow 
IWiegelmann et al.L[2006bh. 

IWiegelmann et al.l d2006bh used a Laplacian-smoothing 
method (first method in Sect. I2.1.2I I in order to reduce the 
amount of small-scale irregularities in the magnetogram. The 
corresponding smoothing functional is calculated as 



i^=2[( A ^ fl -) a+ (^A) a+ ( A ^«) s 



(13) 



A further subfunctional implemented in ppTW measures the 
difference between the preprocessed magnetogram and the orig- 
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inal, observed one, 



L3 = Z K ~~ B( ° bs) ) + K ~ B v° bs) ) + " ^° bs) ) 
p 

(14) 

where B (obs) corresponds to the observed field values. 

The complete minimization functional is then given as the 
sum 

,(TW)_„(TW), ,„(TW). „(TW)r ,„(TW).(TW) 



where fi ( ™\ . . p^ 1 " ' are yet undetermined weighting factors. 
For minimizing this functi onal, ppTW emplo ys a fast and simple 
Newton-Raphson scheme dPress et al.L 1 19891) . 

As an extens io n of the functional given by Eq. dT5l ). 
IWiegelmann et alJ d2008l) implemented another "L5" term, 
which allows to incorporate direct chromospheric observations, 
e.g. Ha images, into the preprocessing scheme. This can im- 
prove the approximation of the chromospheric field, in particu- 
lar the transverse components, from photospheric measurements 
by preprocessing. The ppTW preprocessing scheme has cur- 
rently been implem ented and tested in spherical geometry by 
iTadesse et all d2009l) . 



(TW) 



2.3. Pattern-preserving preprocessing (ppMF) 

In ppMF, the total magnetic force and the total magnetic torque 
are calculated in the same way as in ppTW, namely using Eqs. 
(fTT1 i and <fT~2T >. Additionally, the two quantities are made dimen- 
sionless by normalizing L\ to 



Nr., = 



2 (b^j 



(16) 



values are expected to be small. Also, in the practical extrapo- 
lations an artificial margin where the field vanishes is in general 
added to the magnetogram. 
L^ MF) is normalized to 

NjjMF) = ^ ^]{M„[^ obs) (jt- -n-h,y-n-h),... 

i=x,\.z p 



,Bf bs \x + n ■ h,y + n ■ h)] + B ( ° hs \x,y)} 



(19) 

so that the normalized quantity L^ mf Va^(mf) is ~ 1 for a field of 
maximum roughness. 

The total functional minimized in ppMF is then 



L (MF) = (MF )±L+ (MF)_L L + 

N Lx F2 N L2 H4 



r(MF) 



(MF) ^4 



NjjMF) 



(20) 



with weighting factors /Uj , /4 MF) , an d l^ ¥ \ an d f° r the min- 
imization the me thod of simulated annealing is used (see, e.g., 
iPress et al.L[l989h . 

Most important, differently from ppTW, ppMF does not use 
the subfunctional L3, which provides a global measure of the dif- 
ference between the modified and original photospheric fields. 
Instead, domain borders that must not be overstepped are pre- 
scribed locally at each grid point of the vector magnetogram, 
different for the line-of-sight and transverse field components. 
In this way, modifications induced by the ppMF algorithm into 
the measured field can be consistently limited to an estimation of 
the measurement errors. In the present study (Sect. [3]below) we 
work with overall domain borders of 120 G for the normal field 
and 250 G for the transverse field, but the method may be refined 
by applying domain borders that vary from point to point. 



and normalizing L2 to 



N Ll 



N X Ny 



ZZ v ^W 0h) ) 2 



k=l 1=1 



(17) 



with N x and N r denoting the numbers of grid points in the x 
and y directions (N x ■ N y = N), h the grid spacing (assumed to 

be equal in the x and y directions) and Bj,° bs) the elements of 
the magnetogram matrix. The normalized quantities L\ /Afc, and 
Li/N^ are independent of the units of length and magnetic field 
and both are typically ~ 1 for a field that is neither force-free nor 
torque-free. 

For the smoothing, ppMF uses a windowed-median tech- 
nique (second method in Sect. I2.1.2t . In discretized form, the 
smoothing functional is given as 

4 MF) = Z Tj{ Mn l Bi{x -n-h,y-n-h),... 

i=x,y,z P 

Bi(x + n ■ h,y + n ■ h)] - B t (x, y)f , (18) 

where n is a positive integer number and M„ the median of a rect- 
angular window with (2n + 1 ) ■ (2« + 1 ) grid points centered about 
the point (x,y). The values at the boundary of the magnetogram, 
where the method cannot be applied, are left unchanged; these 



3. Application of the two preprocessing routines to 
AR 10953 

In this section, the two preprocessing routines described in 
Sect. [2] ppTW and ppMF, are applied to a vector magnetogram 
of AR 10953, and the changes in the magnetogram due to the 
preprocessing are examined. After that, the preprocessed mag- 
netograms, as well as the original one, are each extrapolated with 
two differe nt methods, namely the magnetofrictional relaxation 
method of IValori et al J d2007l) and the optimization method of 
IWiegelmann! d2004h . and the resulting fields are analyzed. 

AR 10953 has recently b e en stu di ed using vector mag- 
netograms by lOkamoto et al.1 d2008l) . ISu et alJ d2009l) . and 
iDeRosa et aTT d2009l) . In particular, the main aim of the study 
bv lDeRosa et alJ was to judge the quality of different extrapola- 
tion methods. Here, the attention is focused on the implications 
of using different preprocessing methods. 

3.1. NOAA AR 10953 and data used 

NOAA AR 10953 was a moderately flare-active region with a 
basically bipolar magnetic structure. The vector magnetogram 
employed here is id entical to the one used in the study of 
IDeRosa et al] d2009h . It was derived from the Stokes profiles 
of the two magnetically sensitive Fe I lines at 6301.5 A and 
6302.5 A measured by the Spectro-Po larimeter (SP) instrum ent 
of the Solar Optical Teles cope (SOT: | Tsuneta et all 120081) on- 
board the Hinode satellite dKosugi et al.Ll2007b . The correspond- 
ing scan of the active region began at 22:30 UT on 2007 April 
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30 and took about 30 min. The angular resolution in the East- 
West and North-South directions was 0.32". The 180° ambigu- 
ity in t he transverse field was removed by me ans of the AZAM 
utility dLites et all Il995b iMetcalf et all 120061) . More details on 
the preparation of th e ve ctor magnetogram m ay be found in 
ISchriiver et al.l d2008l) and lDeRosa et all d2009l) . and references 
given in these papers. 

The area of the Hinode/SOT-SP scan covers the central 
part of AR 10953, which was dominated by a leading spot 
of negative polarity. The active region was flare-quiet above 
the C1.0 level until about two days afte r the scan. Images 
from the H i node X-Ray Telescope (XRT; iGolub et all l2007t 
iKano et al 1 l2008h taken around the time of the scan show 
bright loops that probably make visible bundles of coro- 
nal magnetic field lines, see Fig. QI a )- Similar bundles of 
field lines are apparently seen in images from the Extreme 
Ultraviolet Imagers (EUVI), which are parts of the Sun Earth 
Connection Corona l and Heliospheric I nvestigation (SECCHI) 
telescope packages dHoward et alll2008h onboard the two Solar 
Terrestrial Relations Observatory (STEREO) spacecraft, see 
Fig.Hb). 

A significant part of the area of the vector magne- 
togram seems to be magnetically connected to regions out- 
side this area. This is, for instance, also indicated by three- 
dimensional trajectories of loops as stereoscopically determined 
dAschwanden et all l2008h from S TEREO/SECCHI-E UVI im- 
ages (not shown here, see Fig. 1 in lDeRosa et all [2009h . 

In the study of lDeRosa et al.l (l2009h . the vector magnetogram 
was therefore embedded in a larger line-of-sight magnetogram 
with a lower spatial resolution, ob served by the Miche lson 
Doppler Imager (MDI) instrument dScherrer et all 1 19951) on- 
board the Solar and Heliospheric Observatory (SOHO) space- 
craft. The total magnetogram used comprises, in a tangential 
plane to the photosphere, 320 x 320 pixel 2 , with 1 pixel = 
580 km, so that altogether an area of 185.6 Mm x 185.6 Mm is 
covered, whereas the measured full magnetic vector is avail- 
able for a subarea of 100 Mm x 115 Mm. The photospheric 
320 x 320 pixel 2 region in turn is part of an encompassing 
512 x 512 pixel 2 region for which a line-of-sight magnetogram 
was available from the SOHO/SDI measurements, and from 
this 512 x 512 pixel 2 line-of sight magnetogram the potential 
field in a 512 x 512 x 512 pixel 3 cube above the photosphere 
was calculated usin g a Green's function method described in 
IMetcalf et al.l (2008) (with the normal field component set equal 
to zero at the s i de an d top faces of the cube). In the study 
of iDeRosa et alJ (l2009h . the solution domain for the nonlinear 
force-free fields was then a cuboid with a height of 256 pixels 
above the 320 x 320 pixel 2 magnetogram, with the boundary 
conditions at the side and top faces, for the methods that need 
such conditions, given by the potential field calculation for the 
encompassing 512 x 512 x 512 pixel 3 cube, and assuming the 
photospheric field outside the measured vector magnetogram to 
be vertical, i.e., setting the horizontal field components to zero 
there. 

Here, the preprocessing algorithms are applied to the en- 
larged, 320 x 320 pixel 2 vector magnetogram. The subsequent 
analysis of the effects of the preprocessing on the photospheric 
field is only in part done for this enlarged magnetogram and oth- 
erwise restricted to the smaller photospheric region where the 
magnetic vector was measured. 



3.2. Effects of the preprocessing routines on the vector 
magnetogram 

In this subsection, the preprocessed vector magnetograms are 
analysed with respect to the differences they show compared to 
the observed magnetogram and to each other. Namely, we try 
to assess the smoothness of the magnetograms and how well 
they satisfy the conditions of force-freeness and torque-freeness. 
Furthermore, the function a{r) defined by Eq. ([TJ, whose values 
in the level of the photosphere can be calculated from the magne- 
togram, is considered. This function contains much information 
on the topology of the magnetic field lines. In particular, the field 
lines lie in the surfaces a = constant and, thus, a must have equal 
values at any two photospheric points that are connected by a 
field line above the photosphere. Finally, the deviations from the 
original magnetogram produced by the two preprocessing meth- 
ods are compared. 

The minimization functionals, given by Eqs. dT3T ) and d20l >. 
depend on the weighting factors j«' TW) and // MF) , respectively, 
which are free parameters (actually only their ratios count, so 
that for each of the procedures one of the factors can be set 
equal to unity). The optimal choice of the weighting factors 
poses a special problem of the preprocessing. As mentioned in 
Sect. 12.1.21 minimizing the total magnetic force and the total 
magnetic torque on the one hand and the noise level on the other 
hand are competing objectives. If, for instance, the weight of 
smoothing in the minimization functional is set to zero, the sub- 
functionals measuring the forces and torques are generally found 
to be minimized very fast to very small values. If smoothing is 
included, however, the reduction of the forces and torques is im- 
peded. On the other hand, a strong weighting of the forces and 
to rques inhibits the smoo thing (cf. consideration of this problem 
in lFuhrmann et all 120071) . 

For the results presented in the following, the weighting fac- 
tors are chosen such that the competing minimization objectives 
appear properly balanced. The weighting factors for ppTW are 
chosen as ^ = ^ = 1.0, M < TW) = 0.001, m'™" = 0.01, and 
those for ppMF as pf^ = pf F> = ^ MF) = 1.0. ppTW uses a 
normalization of the magnetic field to the mean value of |/? (obs, | 
over the magnetogram and of lengths to the edge length of the 
magnetogram; so the subfunctionals of Z/ TW) and their weighting 
factors are dimensionless but, unlike the corresponding quanti- 
ties in ppMF, not independent of the chosen units of magnetic 
field and length (thus the parameters cannot be compared di- 
rectly). Furthermore, ppMF is carried out with a window size of 
n = 1 for the smoothing and with domain borders of 120 G for 
the normal field component and of 250 G for the transverse field 
components during the minimization. 

Fig. [2] shows outcomes of preprocessing the measured vec- 
tor magnetogram by the two methods. Obviously, ppTW leads 
to a markedly smoother vector magnetogram than ppMF. Small- 
scale irregularities are largely removed by ppTW. The Laplacian 
smoothing used by this method resembles a diffusion process. 
Correspondingly, originally separated maxima or minima may 
be fused, as is seen for the B y component in Fig. [3] where a 
magnification of the central part of the magnetogram region is 
shown. By contrast, though also smoothing the magnetogram, 
ppMF manages to stay close to the observed structures. The 
reason for the different smoothing levels seems to be rooted in 
the different smoothing schemes. The subfunctional L:^^, de- 
fined by Eq. dT8l . decreases very fast during the preprocessing, 
without changing the magnetogram very much. Experiments for 
ppMF with a Laplacian smoothing led to much smoother mag- 
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Fig. 1. Two coaligned images of AR 10953 (with the same 10 degree gridlines drawn on both images for reference), (a) Time- 
averaged and logarithmically scaled Hinode/XKl soft X-ray image and (b) STEREO-A/SECCHI-EUVI 171A image. 



Table 1. Quantitative comparisons of the original and prepro- 
cessed vector magnetograms. a 





Unpreprocessed 


ppTW 


ppMF 


Total magnetic force 


6.71 TO" 2 


2.20 TO -4 


1.48 -10- 6 


Total magnetic torque 


5.50 TO" 2 


9.83 -10- 5 


5.83 -10- 5 


Net magnetic flux M f 


-0.13 


-0.13 


-0.15 


max \a(x, y, z = 0)| 


2.73 


1.22 


2.00 


(\a(x,y,z = 0)\) 


0.057 


0.055 


0.062 



a The total magnetic force, the total magnetic torque, and the net mag- 
netic flux are given in normalized, dimensionless forms (cf. text), and 
a is measured in pixel -1 (1 pixel = 580 km). The total magnetic force 
and the total magnetic torque are calculated for the enlarged, 320 x 320 
pixel 2 magnetogram, the other quantities for the smaller area where the 
magnetic vector was measured. 



netograms, but the decrease rate of the smoothing subfunctional, 



r (TW) 

r (MF) 



was smaller than that of the correponding subfunctional, 
, when applying the windowed-median smoothing. 
A quantitative comparison of the two preprocessing routines 
is given in Table [JJ The values of the total magnetic force and 
and the total magnetic torque shown there were calculated ac- 
cording to Eqs. ( fTTT i and (TT~2T > and subsequently normalized to 
the quantities given by Eqs. dT6b and (fTTl i, respectively. 

ppTW was able to decrease the total force and the total torque 
by two orders of magnitude. The measured total force after ap- 
plication of ppMF is even decreased by four orders of magnitude 
compared to the original total force, and the the measured total 
torque after application of ppMF is roughly by a factor of two 
smaller than that after application of ppTW. 

We argue that the reason for the differences in the capabilities 
to reduce the total magnetic force and the total magnetic torque 
lies in the scheme of the simulated annealing used for the min- 
imization in ppMF, compared to the Newton-Raphson scheme 
employed by ppTW. While ppTW is forced to stop at the nearest 
local minimum it can find, ppMF is able to leave this minimum 
and may find another one. In Fig. |4] which shows the evolutions 
of the total functional L (MF) and of the sum of the normalized 



total force and the normalized total torque in the course of the 
minimization in ppMF, one can see that after roughly 1000 iter- 
ations a local minimum is found, but then this minimum is left to 
find a deeper one after about 2000 iterations. The values of force 
and torque in the first minimum are approximately the ones of 
ppTW. So it seems possible that ppTW was forced to stay at this 
minimum. An alternative possibility is that ppTW was able to 
find the global minimum, but that for the parameters chosen, for 
instance as a result of stronger smoothing, the above minimum 
is not the deepest one. Furthermore, the pathways of the mini- 
mization in the two algorithms will in general be different. 

While for the enlarged, 320 x 320 pixel 2 magnetogram, to 
which the preprocessing is applied, the net flux per unit area 
pixel 2 , 

M ''wikv (21) 

nearly vanishes {unpreprocessed: 2 x 10~ 7 , ppTW: 3 x 10 -7 , 
ppMF: 4 x 10~ 7 ), the vector magnetogram is not flux balanced 
to the same degree, see third row in Table QJ ppTW preserves 
the measured net flux through the vector magnetogram area, 
while ppMF increases it slightly. For the enlarged, 320 x 320 
pixel 2 field of view, both preprocessing routines leave the flux 
balanced. So ppMF moderately redistributes the magnetic flux 
between the area of the vector magnetogram and its surround- 
ings. 

The fourth and fifth rows of Table [TJ and Fig. [5] show how 
the preprocessing routines change the function 



a(x, y, z = 0) 



d x B y {x,y,z = 0) - d y B x {x,y,z = 0) 
B z (x,y,z = 0) ' 



(22) 



Since this quantity depends on all three vector components, and 
is sensitive to noise due to the differentiations needed to cal- 
culate it, it is expected to be particularly strongly affected by 
the preprocessing. The maximum absolute value of a is signif- 
icantly decreased by both procedures (by about 50% by ppTW 
and by about 35% by ppMF, see fourth row in Table [JJ. The 
mean value of |a| over the vector magnetogram area, on the other 
hand, is only weakly changed by both procedures (fifth row in 
Table [JJ. The unpreprocessed a function shown in Fig. [5] (left) 
is highly complex and predominantly negative. ppTW yields a 
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Fig. 2. Original and preprocessed vector magnetograms of AR 10953. From left to right: B x , B y , and B z . From top to bottom: 
Unpreprocessed, ppTW, and ppMF. The magnetic field is measured in G and the length unit is pixel (580 km). 



Table 2. Absolute values of the differences between the prepro- hand, changes the whole magnetogram, but not as strongly as 



cessed and observed fields. 



Preprocessing method ppTW ppMF 



max. (|S, - B xfibs \) [G] 
max. (\B y - B yi0bs \) [G] 
max. (\B Z - B zxlbs \) [G] 
(\B x -B x , obs \) [G] 

(IV 

(IV 



■ Sy.cAs]} [G] 
B : .obs\) [G] 



306.9 
132.8 
215.0 
79.7 
72.9 
93.8 



119.9 
115.7 
116.6 
73.0 
73.4 
81.8 



much smoother a than does ppMF (Fig. [5] middle and right), 
confirming the above observation that ppTW tends to remove 
small-scale structures much more than ppMF. 

Finally, Table [2] shows the absolute values of the differences 
between the preprocessed and observed magnetograms, see also 
Fig. [6] ppTW leads to great changes in strong field areas, but 
leaves weak field areas mostly unchanged. ppMF, on the other 



ppTW. The average changes of B x , B y , and B z over the vector 
magnetogram are nearly identical for ppTW and ppMF. The dif- 
ferences between the two methods may result from the fact that 
in ppTW the deviation from the original magnetogram is con- 
trolled globally via the subfunctional L3, while ppMF uses a lo- 
cal control at each grid point. 



3.3. Effects of the preprocessing on the results of 
extrapolations 

The observed and the two preprocessed maps of the photospheric 
magnetic field vector are now used as input data for extrap ola- 
tions using the optimization method of IWiegelmannl (l2004f) an d 
the magnetofrictional relaxation method of IValori et al.l d2007l) . 

The side and top boundaries are tre ated differently by th e two 
extrapolation methods. The method of IWiegelmannl d2004l) pre- 
scribes the magnetic vector as given by an initial potential field. 
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Fig. 3. By in a magnified cutout of the vector magnetogram of AR10953. From left to right: Unpreprocessed, ppTW, and ppMF. The 
magnetic field is measured in G and the length unit is pixel (580 km). 
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Fig. 5. Contours of the function a(x,y,z = 0). From left to right: Unpreprocessed, ppTW, and ppMF. The length unit is pixel 
(580 km) and a is measured in pixel -1 . 



Our ext rapolations wi t h this method are done, as the extrapola- 
tions in IDeRosa et all d2009l) were done, for a cuboid above the 
enlarged vector magnetogram, and with the boundary condions 
at the side and top boundaries given by the values of the poten- 
tial field calculated for the encompassing 512x512x 5 12 pixel 3 
cube. Contradictions between strong currents in the region above 
the vector magnetogram and the assumption of a potential mag- 
netic field at the side boundaries can be mitigated by placing the 
side boundaries far away from the vector magnetogram area. 



The method of IValori et al.l d2007l) . as used here, applies a 
special kind of open boundary conditions at the side and top 
boundaries, where the normal field component is determined 
from inner field values such as to ensure the solenoidal prop- 
erty of B, and the transverse field is obtained from a fourth-order 
polynomial extrapolation of interior field values to the boundary; 
this precedure works with a buffer layer of three grid points at 
the boundaries. These special boundary conditions allow extrap- 
olations without assuming a potential field at the side and top 
bo undaries. Accordin gly, our extrapolations using the method 
of IValori et al.l (120071) start from the smaller photospheric area 
of the vector magnetogram. That is to say, the solution domain 
for these extrapolations is a cuboid with a height of 256 pix- 
els above the measured vector magnetogram. The necessity of 
modelling the transverse photospheric field outside the vector 
magnetogram area is avoided. Such modelling normally pro- 
duces unphysical current concentrations at and abov e the bound- 
ary of the vector magnetogram. Indeed, using the IValori et"aTI 
method with an embedding in the larger line-of-sight magne- 



togram (IDeRosa et all l2009h leads to worse extrapolation re- 
sults. 

The potential field calculated for the encompassing 512 x 
512 x 512 pixel 3 cube provides the intial field for both types 
of nonlinear force-free field computations. The analysis of the 
extrapolation results is always restricted to a cuboid above the 
area o f the measur ed vector magnetogram, still slightly reduced 
by the IValori et al.l buffer layer. 

For measuring the degree of force -freeness of the extrap- 

c (cf. jVVheatland et al. . 20001: 



ISchriiver et al 1 120061: 1 Valori et all 12007b IDeRosa et all 120091) 



CW si , 



with Jt = (V x B)i, o-i 



\Ji x Bj\ 



(23) 



Isinftl 



where the summation is over the grid points in the considered 
three-dimensional domain V. This is a current-weighted average 
of the sine of the angle 6 between the magnetic field B and the 
current density /, with CW^ n = for an exactly force-free mag- 
netic field. 

Another point of interest is how well the codes fulfill the 
solenoidal condition, and how this is affected by the preprocess- 
ing. Here we use the metric 



Nl 1/1 w M^f \B t \ 
with M being the number of grid points in V. 



(24) 
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Fig. 6. Absolute values of the differences between the preprocessed and observed magnetograms for ppTW (top row) and ppMF 
(bottom row). From left to right: Differences for B x , B y , and B z . The magnetic field is measured in G and the length unit is pixel 
(580 km). 



Similarly, the relative magnetic energy in the volume, de- 
fined as 



" potential,/ 

can be revealing. B po tentiai is the potential field in V whose nor- 
mal component B„ on the boundary dV is identical to that of 
B. For a given distribution of B„ on dV, B po tentiai minimizes the 
magnetic energy content of V, and s > 1 for any non-current-free 
magnetic field B. It has been observed, however, that nonlinear 
force-free extrapolations from observed magnetograms can pro- 
duce fields with s < 1. It has also been noted that this patho- 
logical b ehavior can be corre cted by preprocessing the magne- 
tograms dMetcalf et al.ll2007l) . 

In Table [3] the above metrics are given for the results 
of extrapolations using either the magnetofrictional relaxation 
method or the optimization method and starting from the ob- 
served magnetogram or a magnetogram that was preprocessed 
by ppTW or ppMF. It is seen in the table that for both extrap- 
olation methods, only ppMF improves (i.e., makes smaller) the 
CW s in values of the calculated fields. On the other hand, their 
<|V ■ B\/\B\) values are improved (i.e., brought closer to zero) for 
both extrapolation methods by both ppTW and ppMF. Similarly, 
the relative magnetic energy, e, which is at unphysical values 
below 1 for the fields obtained by the extrapolations starting 
from the observed magnetogram, is raised to values larger than 
1 for both extrapolation methods by both ppTW and ppMF. 
Additionally, we examined the average and maximum values 
of the Lorentz force in the extrapolation volume. ppMF leaves 
these values practically unchanged, while ppTW moderately in- 
creases them. This is largely in agreement with the results for 



the CW s i n value, which stronger weights areas with a high cur- 
rent density. 

Fig. |7] finally, compares current densities for the six extrap- 
olation/preprocessing cases in the form of |V x B\ integrated ver- 
tically over a 12 pixel thick bottom region. For both extrapola- 
tions starting from the unpreprocessed vector magnetogram, the 
obtained current densities show very complex patterns. 

In the case of the magnetofrictional extrapolation method 
(see Fig. [7] top row), the current density obtained when apply- 
ing ppMF appears to be very similar to that for the case with no 
preprocessing. This means that ppMF here leads to a solution 
close to that obtained without preprocessing, but with a physical 
magnetic energy content and in better agreement with the force- 
free and solenoidal conditions (cf. Table 0. The field calculated 
by magnetofrictional extrapolation after applying ppTW, on the 
other hand, seems to be a much smoother and less structured 
version of the one calculated without preprocessing, at the same 
time being more physical in having an energy content in excess 
of the potential field value and in better fulfilling the solenoidal 
condition. 

For the extrapolations using the optimization method (see 
Fig. [7] bottom row), the two preprocessing algorithms seem 
to have broadly similar effects on the current density. If one 
compares the cases with no preprocessing and with ppTW ap- 
plied, one can again see that ppTW smoothes the current density 
strongly. In the case of preprocessing with ppMF the smooth- 
ing effect is also observable, though to a lesser degree. This 
may have to do with the fact that, independently of the kind 
of preprocessing applied, the current density distributions ob- 
tained by extrapolations with the optimization method show 
weaker contrasts than those obtained by magnetofrictional ex- 
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Table 3. Comparison of extrapolation results for different preprocessing and extrapolation methods. 3 



Extrapolation 


Boundary map 


cw sm 


(\j x 0|) 


max(L/x0|) 


(|V-0|/|0|) 


s 


Magneto- 


unpreprocessed 


0.12 


0.12 


62.1 


10.8 


0.67 


frictional 


ppTW 


0.15 


0.21 


88.3 


6.4 


1.14 


method 


ppMF 


0.08 


0.14 


62.2 


3.6 


1.12 


Optimi- 


unpreprocessed 


0.35 


0.23 


97.6 


7.2 


0.87 


zation 


ppTW 


0.44 


0.30 


134.3 


3.8 


1.05 


method 


ppMF 


0.29 


0.20 


98.9 


1.2 


1.02 



a CW s i„ and e are dimensionless, (|/ x 0|) and max (jj x 0|) are measured in 10 3 G 2 pixel 
580 km). 



and (|V • 0|/|0|) is measured in 10 7 pixel ' (1 pixel 
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Fig. 7. Absolute value of current density, |VxZJ|, integrated vertically over a 1 2 pixel thick bottom layer, in G. Top row: Extrapolation 
with magnetofrictional method. Bottom row: Extrapolation with optimization method. From left to right: Unpreprocessed, ppTW, 
and ppMF. The length unit is pixel (580 km). 



trapolations. Comparing the figures in Table [3] from the point 
of view of the solutions' consistency, we can see that the mag- 
netofrictional extrapolations are closer to force-free (by a fac- 
tor of three in CWsin), while the optimization method led to 
more solenoidal reconstructions. The optimization method is rel- 
atively little affected by the energy in the small scales of the un- 
preprocessed magnetogram, while the magnetofrictional method 
attained higher free energies with preprocessed magnetograms. 
(We do not extend further these observations since the compari- 
son of different extrapolation techniques is not our subject here). 

4. Conclusions 

We have compared and discussed the two so far existing meth- 
ods for t he preprocessing of solar ve ctor magnet ograms, namely, 
those o f TWiegelmann et al.l (l2006bl ppTW) and iFuhrmann et all 
(120071 ppMF). These methods make the magnetograms more 
suitable for nonlinear force-free extrapolations into three- 
dimensional magnetic fields in the chromosphere and corona. 



Both methods follow the same strategy, namely, to minimize a 
functional L of the photospheric field values such as to simulta- 
neously make small the total magnetic force and the total mag- 
netic torque on the volume considered and the amount of small- 
scale noise in the photospheric boundary data. 

The two methods differ in their ways of minimizing the 
functional L. While ppTW employs a Newton-Raphson scheme, 
ppMF uses the method of simulated annealing. Furthermore, the 
task of reducing the small-scale noise, or smoothing, is solved 
differently. ppTW applies a minimization of the absolute value 
of AB(x, y,z- 0), which corresponds to algebraically averaging 
the photospheric field, while in ppMF the smoothing is reached 
by a windowed-median averaging. Finally, the degree of devia- 
tion of the modified fields from the original photospheric field 
during the minimization is controlled in different ways by the 
two methods. In ppTW a special subfunctional that measures the 
deviation is included as an additional minimization target in L. 
This method of control is integral in character. ppMF, on the 
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Fig. 4. Evolution of the total functional L (MF) as given by 



Eq. <H)J (top) and of the sum p.^ Li / N L , + pf ¥) L 2 IN Ll of the 
normalized total force and the normalized total torque (bottom) 
during the minimization in ppMF on a logarithmic scale. Both 
quantities are calculated for the enlarged, 320 x 320 pixel 2 mag- 
netogram. 



other hand, uses a local control at each grid point in the form of 
an interval within which the field values have to stay. 

The two preprocessing methods were applied to a vector 
magnetogram of the recently observed active region NOAA AR 
10953, which had already been the target in a comparative study 
of different extrapolati on methods for nonl inear force-free mag- 
netic fields reported in lDeRosa et al.l d2009l) . Both preprocessing 
methods managed to significantly decrease the magnetic forces 
and magnetic torques that act through the magnetogram area and 
that can cause incompatibilities with the assumption of force- 
freeness in the solution domain. The force and torque decrease 
was stronger for ppMF than for ppTW. 

Both methods also reduced the amount of small-scale irreg- 
ularities in the observed photospheric field, where ppTW led to 
a markedly smoother magnetogram than ppMF, in accordance 
with the aim of ppTW to mimic the expansion of the solar mag- 
netic field between photosphere and chromosphere. The aver- 
age deviations of the preprocessed magnetograms from the ob- 
served magnetogram are nearly identical for ppTW and ppMF, 
where ppTW led to greater changes in strong-field areas, leav- 
ing weak-field areas mostly unchanged, while ppMF weakly 
changed the whole magnetogram, thereby better preserving pat- 
terns present in the original magnetogram. Similarly, for the 



function a(x,y, z = 0) ppTW yielded a much smoother distri- 
bution than ppMF. 

The original magnetogram and the two preprocessed mag- 
netograms were used as input data for nonlinear force-free field 
extrapolations by means of two differe nt methods, namely the 
magnetofrictional relaxation method of Valori et all d2007l) and 



the optimization method of IWiegelmannl ( 20041) . Both ppTW 
and ppMF corrected a pathological property of the fields cal- 
culated by the extrapolations from the observed magnetogram, 
namely they raised the magnetic energy content of the extrapo- 
lated fields from values below to values above that for the po- 
tential field with the same normal component on the boundary. 
Also, the fields calculated from the preprocessed magnetograms 
fulfill the solenoidal condition for B better than those calculated 
without preprocessing. Finally, some effects of the preprocess- 
ing, as those on the current density in a bottom layer, seemed to 
be influenced by the the specifics of the employed extrapolation 
method. 
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